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Summary 


For  a  sample  of  individuals  from  an  animal  or  human  population 
under  observation  in  a  clinical  trial  or  life  test,  mortality  rates  are 
defined  for  each  risk  to  which  the  population  is  exposed.  By  considering 
only  the  crude  risk  functions  and  cause-specific  failure  rates  in  a  basic 
competing  risks  model,  these  mortality  rates  are  shown  to  have  an  asymptotic 
normal  distribution.  An  expression  for  the  asymptotic  correlation  be¬ 
tween  a  pair  of  mortality  rates  is  thus  obtained  and  a  necessary  and 
sufficient  condition  for  their  asymptotic  independence  is  investigated  in 
some  general  situations  with  reference  to  only  the  cause-specific  failure 


rates  and  the  survival  function  of  the  individuals. 
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1.  Introduction 


•  1 

We  shall  consider  here  a  basic  statistical  model  in  which  an  animal 
or  human  population  is  under  insult  by  a  totality  of  k  mutually  exclusive 
competing  risks.  We  are  concerned  with  experiments  in  which  a  subpopulation 
of  n  units  is  under  continual  surveillance  from  the  onset  with  the  time 
and  type  of  occurrence  of  the  responses  monitored,  whenever  possible,  for 
each  unit.  The  term  'response'  is  given  a  wide  interpretation  and  encompasses 
a  broad  class  of  response  data  from  clinical  and  epidemiologic  studies  in 
which  the  nature  of  the  observed  response  can  be  attributed  to  exactly  one 
of  several  causes.  In  a  variety  of  such  investigations  a  group  of  indivi¬ 
duals  may  be  followed  with  interest  in  a  specific  cause  of  the  observed 
response.  For  instance  in  a  laboratory  animal  experiment  concerned  with  the 
exposure  to  certain  carcinogenic  pollutants  the  terminal  response  may  be 
the  appearance  of  a  tumor  or  death  caused  by  one  of  several  organotropic 
carcinomas  --  intestine,  bladder,  lung,  liver  —  but  concern  may  be  focused 
on  only  one  specific  carcinoma.  Experiments  of  this  type  are  often  com¬ 
plicated  by  the  presence  or  need  for  censoring  designs.  Limitations  on  time 
or  economic  considerations  may  prevent  the  luxury  of  prolonged  experimenta¬ 
tion  beyond  specified  constraints  or,  in  some  instances,  a  desire  for  a 
preliminary  analysis  of  the  current  accumulated  information  in  a  study 
necessarily  entails  incomplete  data.  On  the  other  hand  withdrawals  from  a 
clinical  study  and  subsequent  loss  to  follow-up  must  be  considered  in  some 
analyses. 

In  this  paper  we  consider  the  mortality  or  response  rate  of  the  com¬ 
peting  factors  evaluated  from  data  recorded  in  a  clinical  trial  or  life 
test  in  which  n  specimens  are  observed  either  throughout  their  lifetime 
or  over  a  duration  of  time  (Type  I  censoring).  The  statistical  model  which 
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we  describe  will  also  permit  for  designs  in  which  items  are  withdrawn  from 
observation  at  random  (Random  censoring). 

The  mortality  rate  for  a  risk  type  is  defined  as  the  ratio  of  the 
observed  frequency  of  failures  among  the  n  specimens  attributed  to  that 
risk  to  the  total  observation  time  for  all  units.  However,  as  the  occurrence 
of  a  failure  by  one  cause  necessarily  precludes  failare  by  the  other  risks 
these  mortality  rates  will  in  general  be  dependent  and  thus  any  statistical 
procedures  based  on  them  for  inferences  on  a  particular  risk  type  should 
take  into  account  the  extent  of  the  possible.influence  of  the  remaining 
risks.  We  shall  show  that  under  very  general  conditions  the  vector  of 
mortality  rates  for  the  k  competing  risks  operating  on  a  random  sample  of 
n  specimens  has  a  limiting  k-variate  normal  distribution-.  An  explicit 
expression  for  the  asymptotic  correlation  of  a  pair  of  mortality  rates  is 
obtained  and  the  sign  of  this  correlation  is  seen  to  be  related  to  the 
monotonicity  of  the  cause-specific  failure  rates  of  the  individual  risks 
(when  all  risks  are  operative).  Examples  are  cited  where  this  correlation 
is  zero  thereby  ensuring  asymptotic  independence  of  these  two  rates.  We 
also  explore  this  in  the  proportional  hazards  model  of  Chiang  (1968). 

The  evaluation  of  mortality  rates  and  their  use  in  some  survival 
analyses  have  been  studied  by  several  researchers  in  recent  times.  Tolley 
and  Norman  (1979),  Vaeth  (1979)  and  Zelen  and  Dannemiller  (1961)  draw 
attention  to  the  possible  influence  of  competing  risks  on  a  specific 
mortality  rate  under  study.  Our  results  based  on  very  general  assumptions 
(see  Section  2)  indicate,  at  least  asymptotically,  the  extent  of  this 
dependency  and  now  allow  for  procedures  under  broader  conditions.  Byar 
(1977)  has  advocated  the  use  of  mortality  rates  evaluated  from  various 
categories  defined  by  prognostic  variables  in  some  medical  statistical 
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investigations  and  has  given  testing  procedures  for  the  homogenity  of  these 
rates  on  the  assumption  that  the  response  times  of  the  competing  variables 
follow  exponential  distributions.  A  more  general  picture  now  emerges  from 
our  results. 


2.  Preliminaries 

We  consider  a  random  sample  of  n  specimens  under  observation  from 
the  onset  until  time  t.  For  each  subject  we  record  either  the  time  of 
failure  (death)  Y  and  the  associated  cause,  or  the  censoring  time  t.  The 
variable  J  €  (l,...,k)  labels  the  specific  cause  of  failure,  if  observed 
before  censoring.  For  the  purpose  of  our  analysis  it  is  not  necessary  to 
consider  latent  (conceptual )  life  times  corresponding  to  the  various  risks 
and  a  joint  distribution  of  these  lifetimes  as  is  customary  in  seme  situa¬ 
tions  (for  example,  Moeschberger  and  David  (1971),  Hoel  (1972),  Gail  (1975) 
and  Elandt-Johnson  (1976)).  Tsiatis  (1975)  and  Peterson  (1976)  have  dis¬ 
cussed  the  serious  identifiability  problems  that  rise  in  connection  with 
assumptions  about  this  joint  distribution  of  conceptual  life  times  when 
employing  data  of  the  type  (Y,J).  Accordingly  we  shall  only  deal  with  the 
absolutely  continuous  functions  Q. ,  called  crude  risks,  defined  by 

Qi ( t )  =  PCY  <  t,  J  =  13,  t  >  0,  i  =  1 . k  (1) 

and  the  cause-specific  failure  rates  g^  given  by 

g. (t)  =  lim  P[t  <  Y  <  t  +  At,  J  =  i|Y  >  t]/At,  i  =  l,...,k.  (2) 

1  At-*0 

We  may  interpret  g^(t)  as  the  instantaneous  probability  of  failure  by 
cause  i  at  time  t  when  all  risks  are  operative,  given  survival  from 
all  risks  up  to  time  t.  In  some  contexts  the  g^  have  also  been  referred 
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to  as  the  rates  of  decrement,  hazard  rates, forces  of  mortality  and  intensity 
rates.  Now,  by  virtue  of  (1)  and  (2)  we  have 


where 


Qi (t: )  =  j  g^yJSCyJdy,  t  >  0 


(3) 


S(t)  =  P[Y  >  t],  t  >  0  (4) 

is  called  the  survival  function  of  the  units.  Since  the  k  risks  are 
mutually  exclusive  and  exhaustive  one  has 

k  k 

l  Qi(t)  =  1  -  S(t)  and  l  Q-H  =  1  .  (5) 

1=1  i=l  1 

2 

Throughout  this  paper  we  make  the  mild  assumption  that  EY  >  0  and.  EY  <  <*>. 

Suppose  the  n  subjects  are  monitored  over  the  time  interval 
[0,t).  The  total  observation  time  is 

Tn(t)  =  l  Y  a  t  (6) 

n  £=T  * 

and  the  number  of  recorded  failures  due  to  cause  i  is 

°n,i(t>'  j,  <7) 

Here  x  a  y  =  min(x,y)  and  1(A)  denotes  the  indicator  of  a  set  A.  In 
some  contexts  Tn(t)  has  been  designated  the  'time  on  trial’  or  'person 
years  at  risk'  of  the  sample.  The  mortality  rate  for  the  cause  i  evaluated 
from  this  information  recorded  in  the  interval  [0,t)  is  defined  as 

r,M(t>  ■  °n,1(t>/Vt>-  <8> 

We  are  interested  in  the  vector  of  mortality  rates  designated  by 

-n^  =  ^rn,l^t^’*"*rn,k^t^’ 


(9) 
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If  all  units  are  monitored  over  their  entire  life  time  the  corresponding 
mortality  rates  can  be  obtained  from  the  above  by  setting  t  equal  to 
infinity.  To  deal  with  random  withdrawals  we  suppose  each  individual  has 
its  own  withdrawal  time  Z  and  so  corresponding  to  (6)  and  (7)  we  may  de¬ 
fine 


*  5 

Tn  =  I 
£=1 


vz 


£•  Dn.i"  j,  I[Y^V  Hi<k- 


(10) 


In  this  formulation  we  may  consider  withdrawals  as  a  competing  risk  with 
its  own  survival  time  Z.  With  this  interpretation  the  situation  reduces 
to  that  of  observing  the  sample  until  all  units  have  responded. 


3 •  pistrjbutjpn. of. mortal  ity,  rates,  and  discussion 
In  view  of  our  definitions  (1)  and  (4)  we  obtain  from  (6)  and  (7), 
for  each  i,  1  <.  i  <.  k  and  t  0 

E  DnJ(t)  =  £  P[Y£  <  t,  J£  =  i]  =  n  Qi(t),  (11) 

E  T  (t)  =  l  f1  P[Y  >  yidy  =  n  f  S(y)dy.  (12) 

n  £=1  Jo  1  J0 

Therefore  using  (3)  we  define  the  vector  ^(t)  =  (y-j  (t),. ..  .YjJt))  where 


m  E°n.i(t)  J0 

Yi(t)  ~  E  T  (t)  ‘ 


9i(y)S(y)dy 


,  1  <  i  <  k,  t  >  0  . 


(13) 


S(y)dy 


We  shall  also  deal  with  a k  *  k  matrix  £( t)  whose  ( i , j ) th  element  o..(t) 

rst  I  J 

is  given  by 

°ij(t)  =  {f  6ij  -  cj(t))  +  Qj(t)(«s  6^  -  C.(t))},  (14) 


1  <  i,  j  <  k,  where  6...  is  the  Kronecker  delta  (6^.  =  1,  all  i  and 


6^  =  0  otherwise)  and  C •  ( t )  is  defined  by 

•  J 


- 

ft 

y  g^yNyJdy  ( 
o  1  J 

9,-  (y)S(y)dy)( 
0  1 

yS(y)dy)  -< 
0 

- 

s(y)dy  .  ( 

0  j 

S(y)dy)2 

0 

- 

Proposition .  With  t  held  fixed  and  y(t),  £(t)  defined  by  (13)  and  (14) 
respectively 

LCnls(£n(t)  -  x(t))]  Nk(ft,  £(t))  as  n  ■>  °°  .  (16) 

Corollary.  With  the  interpretation  made  on  random  withdrawals 

LCnJs(r(~)  -  x(°°))3  -*•  \(g»  £(“))  as  n  -*•  “• 

The  proof  of  the  proposition  will  be  taken  up  in  the  next  section. 

We  devote  the  remainder  of  this  section  to  some  remarks  and  examples. 

Remark  1 .  From  (14)  we  have  that  the  asymptotic  correlation  between  the 
mortality  rates  r  (t)  and  r - (t)  is 

*  J 

-  (Qi(t)CJ.(t)  +  Qj(t)Ci(t)) 

/Qi'(t7QJTt)(l  -  2C.(t))(l"-  ?C'j(t)) 

and  thus  (16)  provides  a  necessary  and  sufficient  condition  for  the  asymptoti 
independence  of  these  rates,  namely 

Qi(t)Cj(t)  +  QJ(t)C1(t)  =  0. 

An  analogous  statement  can  be  made  in  the  case  where  experimentation  is 
monitored  throughout  the  life  time  of  all  units  of  the  sample. 

RemaHc_2.  We  shall  show  that  the  C-(t)  of  (15)  may  be  interpreted  as 
covariances.  To  see  this  let  us  define  a  nonnegative  random  variable  Wt 
whose  probability  density  function  is  S(x)/(/J  S(y)dy),  for  x  €  [0,t] 
and  zero  otherwise.  Then  (15)  may  be  rewritten  as  C.(t)  =  Cov(Wt,gi (Wt) ). 
Hence  with  t  held  fixed,  Ci (t)  >  0  (respectively  <  0)  according  as  the 
cause-specific  failure  rate  x  -*  g.j(x)  is  nondecreasing  (respectively 
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nonincreasing)  in  T0,t).  In  particular  C^t)  =  0  if  x  -*■  g ^  ( x )  is  a 

constant  function  on  [0,t).  Note  that  C^t)  is  never  zero  if  g^  is 

(strictly)  increasing  or  decreasing  over  [0,t).  Similar  remarks  hold  in 

the  situation  in  which  the  sample  is  monitored  until  all  units  have  failed. 
Remark  3.  Now  if  for  a  pair  of  risks  i,j  the  failure  rates  g^, 

are  both  nondecreasing  (respectively  nonincreasing)  in  [0,t),  the  asymptoti 

correlation  of  the  mortality  rates  I\(t),  F ^ ( t )  will  be  nonpositive 

(respectively  nonnegative).  In  particular  when  both  g.,  g.  are  constant 

*  J 

functions  on  [0,t),  this  correlation  is  zero  and  hence  from  (16)  this 

ensures  the  asymptotic  independence  of  the  rates  (t)  and  fj(t).  It  may 

also  be  noted  that  if  all  g-  are  constant  over  [0,°°),  then  the  hazard 

1  k 

rate  corresponding  to  the  survival  time  Y  is  also  constant  (=  l  g.), 

1=1  1 

whence  Y  must  have  the  exponential  distribution. 

We  present  an  example  in  which  the  mortality  rates  are  asymptotically 
independent  even  though  the  corresponding  failure  rates  are  nonconstant. 
Example  1.  Consider  a  life  test  in  which  n  items  are  observed  throughout 
their  life  time.  We  specify  a  competing  risks  model  with  crude  risk  func¬ 
tions  given  by 


Q.j (t)  =  X.(l  -  exp(-Ait)),  t  >  0,  i  =  1,2, 


with  X-j  >  X2  >  0  and  X-|  +  *2  =  ^  *  Then  ^rom  ^  anc*  we  have 

S(t)  =  X1  exp(-X-|t)  +  X£  exp(-X2t),  t  ^  0 
g^(t)  =  exp(-X1-t)/S(t),  t  >  0,  i  =  1,2. 


From  the  Corollary  it  will  follow  that  the  mortality  rates  r^00),  I^f00) 
are  asymptotically  independent  provided  (Q-j  (“jCgf")  +  (°°) )  =  0. 

From  (15)  we  have  in  general 


Cj  («) 


E[YI[J  =  i  ]] 
EY 


^  Qj(°°) 


For  our  example  a  straightforward  calculation  shows 
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EY  =  S(y)dy  =  2; 

J0 

EY2  =  2  f  y  S(y)dy  =  2U"1  +  A"1). 

r°° 

Also  Q.(°°)  =  A-  and  E(YI(J  =  i))  =  y  Q!(y)dy  =  1;  i  =  1,2.  Hence 
1  JO 

CtH  =  Ml  -  x^1), 

c2H  =  Mi  -  XgA”1 ) 


and  therefore  (Q-j  (“OCgf00)  +  (^(“K]  (°°) )  =  0.  We  note  that  g^  is  strictly 
decreasing  and  g^  strictly  increasing  whence  in  view  of  Remark  2, 

C^(»)  <  0  and  C2(°°)  >  0  as  is  the  case. 

Remark  4.  From  our  assumptions  the  total  failure  (hazard)  rate  corresponding 

k 

to  the  survival  time  Y  is  h(t)  =  l  g-(t)  with  h(t)  =  -S(t)/S'(t). 

i=l  1 

Now  suppose  Chiang's  proportionality  assumption  holds,  namely  g.(t)  =  c-h(t), 

k  1  1 

where  the  c-  (>  0)  satisfy  £  c.  =  1.  Therefore  if  Y  has  increasing 
1  i=l  1 

failure  rate  (IFR)  then  h  and  so  all  gi  will  be  increasing  and  so 
C.(t)  >  0  for  all  i  and  t.  Hence  any  pair  of  mortality  rates  will  be 
negatively  correlated.  Likewise  if  Y  has  decreasing  failure  rate  (DFR) 
all  pairs  of  mortality  rates  will  be  positively  correlated.  This  holds  for 
any  truncation  time  t. 

Again  from  (15),  we  obtain  after  some  simplification 


C  =  l  C.H  =  (1 
i=l  1 


Under  the  proportional ity  assumption  =  c^C,  for  all  i.  Thus  if  Y 

has  increasing  failure  rate  average  (I FRA)  then  C  >  0  and  again  each  pair 
of  mortality  rates  r.(°°),  r.(°°)  will  be  negatively  correlated.  Of  course 

1  J 

if  Y  has  IFR  it  has  necessarily  IFRA.  Similar  statements  can  be  made  when 
Y  has  decreasing  failure  rate  average  (DFRA). 
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Example  2.  Suppose  the  crude  risks  are  specified  by 

Q.-(t)  =  f1  ya']e'y/%,  t  >  0,  i  =  l,...,k 

1  J0 

k 

where  a,8>  all  pH  >  0  and  without  loss  of  generality  we  take  [  p.  =  1 

i=l  1 

Here  r  denotes  the  gamma  function.  The  survival  time  Y  has  the  gamma 
distribution  with  parameters  a, 13.  We  get 


S(t)  =  {Bar(a)}-1  f  ya_1e  y/2dy 


-CO 

g^t)  =  p.  (1  +  t"1y)a'1e"y/edy  =  p.h(t),  t  >  0 
'  *  1 


where,  as  before,  h  is  the  failure  rate  of  Y.  Thus  the  proportionality 
assumption  holds.  Furthermore  Y  has  IFR  or  DFR  according  as  a  >  1  or 
a  _<  1 .  For  a  =  1,  h  is  constant.  A  simple  calculation  will  show 
C-j(°°)  =  4p^(1  -  a""*)  and  so  unless  a  =  1 ,  no  two  mortality  rates  will  be 
asymptotical ly  independent. 


4-  Ecfiaf^Qf^EmBOsition 

From  (6),  (8)  and  (13)  one  obtains  for  each  1  £  i  <  k,  t  _>  0 

n!S(rn,i(t)  -  Yj (t) )  »  {n"1Tn(t)}"1{n"i5  l  Zu}  (17) 

1 

where  Z =  I C Y^  <  t,  =  i]  -  Y-jUHY^  a  t}.  Now,  for  each  i,  1  <  i  _<  k, 
the  Z^,  1  <_  l  _<  n  are  independent  and  identically  distributed  (iid)  random 
variables.  Also  EfZ.^)  =  0.  We  need  to  compute  Cov(Z . ^,Zj^) .  For 
simplicity  in  script  write  =  IFY^  <  t,  =  i].  Then 


zUzij  *  -  Vt)(vi  A  t)nj*  ■  ',j(t)fV  A  t,nu  +  ViUlYjUHVj  a  t)‘ 

08) 
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But  Wil =  °’ if  1  *  j 

=  ni£  if  i  =  j. 


Also  E{ (Y£  A  t)ni;}  =  j0  y  Qi (y )dy ,  1  <  i  <  k  . 

Hence  from  (18)  through  (20),  if  i  f  j, 

c°v(zi£,zj£)  =  y  Qj (y)dy)  -  Yj(t ){f*  y  Qj(t)dy) 


+  2Yi(t)Yj(t){J  y  S(y)dy}  .  (1 

Recall  again  (13)  and  (15).  Then  after  some  routine  manipulations  (21) 
can  be  rewritten 


Cov(zu,zj£)  =  -(Q.j  (t)Cj(t)  +  Qj(t)Ci(t)}  (e 

for  all  t  >  0  and  1  <  i,  j  <  k,  with  i  f  j.  But  from  (18)  and  (19) 


Var^Z>£^  E^nit  “  2"Yi ( t) ( Y&  a  t)ni£  +  Yi(t)(Y£  a  t)  } 


Cov(Zu,Zu) 


=  Qi (t)  -  2yi (t)  y  Ql(y)dy  +  ZY^tjj  y  S(y)dy 
*  Q i ( t ) ( 1  -  2C-(t)),  1  <  i  <  k,  t  >  0  .  (23) 

To  establish  the  proposition  consider  the  random  variable 

k  n 

*n(t)  =  E  Vn%(rn,i(t)  ~  Y-j (t) ) }  =  {n"1Tn(t)}‘1n^  £  Z*  (24) 

★ 

where  the  a^'s  are  arbitrary  constants  and  the  sequence  (Z£:  1  <_  i  <  n} 
is  given  by 

*  i 

Z4  =  i^1  aiZi i  ' 

Hence  by  the  Central  Limit  Theorem  n~2  7  Z  ■*  U(0,od)  with 

4*1  £  v 
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2 

a  =  l  l  a. a-  Cov(Z.0,Z..).  Also  by  the  Strong  Law  of  Large 
i=l  j=l  1  l 

i  rt 

S(y)dy.  Thus  (24)  leads  to  the  convergence 


Numbers . 


Vs 


X„(t)  N(0,( 


S(y)dy)"2o2). 


(25) 


Writing  a  =  (a^,...,a(c)  and  recalling  (14)  we  have  shown  in  (25)  that 
£{n^(£n(t)  -  x(t)))aT  ->  N(0,a  l  aT).  So  since  a  is  an  arbitrary  vector, 
the  proposition  is  proven.  The  Corollary  can  be  established  along  exactly 
the  same  lines. 


5.  Concluding  Remarks 

In  this  paper  we  have  dealt  with  mortality  rates  evaluated  in  a  life 

test  or  clinical  trial  under  a  basic  competing  risks  model.  The  asymptotic 

.  multivariate  nQrmal  distribution  that  we  have  obtained  for  the  vector  of 

these  mortality  rates  also  provides  for  inference  on  the  parameters  of  our 

model.  For  instance  confidence  sets  for  the  y(t)  can  be  supplied  and 

tests  of  hypotheses  on  ^(t)  obtained.  We  may  also  consider  the  stochastic 

process  t  ->  n2(£n(t)  -  y(t)),  t  >  0  as  an  element  in  the  function  space 

D  CO,00)  (the  k-fold  cartesian  product  of  D[0,»)  space  of  right  continuous 

real  functions  on  [0,°°)  with  left  hand  limits)  equipped  with  a  suitable 

topology,  and  thus  an  invariance  principle  derived  paralleling  our  Proposition. 

Furthermore  appropriate  stopping  mechanisms  defined  in  terms  of  £n(t)  can 

be  considered  as  for  example  the  rule  t  =  inf{t  >  0:  min  F  . (t)  >  a  }, 

l<i_<n  11,1  n 

where  the  an  are  some  constants. 
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